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ABSTRACT 

The results of a spectroscopic programme aiming to investigate the kinematics of the local 
interstellar medium components towards the Southern Coalsack and Chamaeleon-Musca dark 
clouds are presented. The analysis is based upon high-resolution (R k, 60 000) spectra of 
the insterstellar Nal D absorption lines towards 63 B-type stars {d ^ 500 pc) selected to 
cover these clouds and the connecting area defined by the Galactic coordinates: 308° ^ I ^ 
294° and —22° ^ & ^ 5°. The radial velocities, column densities, velocity dispersions, 
colour excess and photometric distances to the stars are used to understand the kinematics and 
distribution of the interstellar cloud components. The analysis indicates that the interstellar gas 
is distributed in two extended sheet-like structures permeating the whole area, one at d ^ 60 
pc and another around 120-150 pc from the Sun. The nearby feature is approaching to the 
local standard of rest with average radial velocity of —7 kms^^, has low average column 
density logiVNai~ 11.2 cm~^ and velocity dispersion b « 5 kms^^. The more distant 
feature has column densities between 12.3 ^ log N-^^ i ^ 13.2, average velocity dispersion b 
« 2.5 kms^^ and seems associated to the dust sheet observed towards the Coalsack, Musca 
and Chamaeleon direction. Its velocity is centered around kms^^, but there is a trend for 
increasing from —3 kms^^ near 6 = 1° to 3 kms^^ near 6 ~ —18°. The nearby low column 
density feature indicates a general outflow from the Sco-Cen association, in aggreement with 
several independent lines of data in the general searched direction. The dust and gas feature 
around 120 - 150 pc seem to be part of an extended large scale feature of similar kinematic 
properties, supposedly identified with the interaction zone of the Local and Loop I bubbles. 
Assuming that the interface and the ring-like volume of dense neutral matter that would have 
been formed during the collision of the two bubbles have similar properties, our results rather 
suggest that the interaction zone between the bubbles is twisted and folded. 

Key words: Stars: distances - ISM: clouds - ISM: individual objects: Southern Coalsack - 
Chamaeleon-Musca - Loop I - Local Bubble 



1 INTRODUCTION 

The present work is part of a comprehensive investigation of the 
local interstellar medium (ISM) towards the Southern Coalsack and 
the Chamaeleon-Musca dark clouds. A comparison of the colour 
excess E{b — y) vs. distance diagrams for the Chamaeleon-Musca 
complex (Franco 1991) and the Southern Coalsack (Franco 1989) 
shows remarkable similarities. The jump of the colour excess to 
higher values occurs approximately at the same distance, and the 
observed minimum value of this rise is almost the same: AE{b — 
y) ~ ff!^ 100. Although the clouds are apart by more than 15° these 
facts suggested that they might be dense condensations embedded 
in an extended interstellar structure. 

In order to investigate the possible physical association be- 
tween these two clouds, a new photometric investigation was car- 
ried out. The campaign resulted on very precise uvbyP data for 



1017 stars covering these clou ds and the connecting are a: 307° ^ 
I ^ 294° and -20° ^ 6 =^ 5° JCorradi & Francoll995h . 

Analysis of the colour excess vs. distance diagrams for these 
stars indicated that exists a foreground region with very little red- 
dening {E{b — y) 0'"006), that is bound, at the distance range 
of 150 ± 30 pc, by a transition region where the colour excesses 
suddenly increase to a mininum value of E{b — y) ~ CP^050 for the 
most diffuse lines of sight and E{b — y) m 0^100 for the densest 
ones. Beyond this transition region the diagrams suggest the ex- 
istence of a second low column density volume, as the complete 
range of colour excess from ~ 0^050 to 0™275, shown by the stars 
in the very narrow distance slot centered on 150 pc, remains un- 
changed for at least another 350 pc (Corradi et al. 1997, hereafter 
Paper I). 

The onset of this minimum colour excess for the whole sur- 
veyed area suggests that the obscuring material is distributed in an 
extended interstellar sheet-like structure. The existence of such an 
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absorbing feature, at a distance identical to the molecular clouds, 
also suggests that Coalsack, Chamaeleon and Musca can be dense 
condensations embedded on the diffuse medium composing the 
sheet-like feature. 

It has been also noted in Paper I that the jump in E(h — y), 
caused by the dust feature, seems to have a dependence on the 
galactic latitude that can be represented by [E{b — 2/)min,6] = 
[0'"050;0°] [0'"100;-8°] [ff"150; -15°]. The quoted 
dependence might indicate, if we postulate that the dust sheet is 
roughly perpendicular to the galactic plane, that the sheet does not 
have the same optical thickness, because E{b — i/)min vs. (6) does 
not follow a simple E{b — y)sec(b) law. On the other hand, it 
might also be the result of approaching the tangential point of a 
warped sheet-like structure with the same column density, curved 
away from the Sun. 

As pointed out in Paper I, when viewed in connection to 
the other data on the local ISM, the existence of these two low- 
reddening volumes has led to the idea that the dust sheet could be 
part of a large scale structure, probably related to the interface of 
the Local and Loop I Bubbles. Such an interface of dense mat- 
ter may have resulted from the compression of the ISM by the 
action of energetic events of the Sco- Cen OB association from 
the far side (Weaver 1979; Iwan 1980: Crawford 1 1 99 it Ide GeusI 
Il992i) and either by a supernova explosion near the Sun or by 
what ever created the local low density region from the near side 
(^e.g. 'Cowie & Soneaila 1986; Cox & Revnolds 1987; Bochkarev* 

[1987; Frisch & York 1986; Gehrels & Chen 1993; Hartauist 1994; 
Bruhweiler 1996). If the dust sheet is part of a small bubble, it could 
be associate d to the LCC-shell pr oposed bv yde Geusi (.1993) . Note 
however that lGenova et alJ il997h find that gas is moving out from 
Sco-Cen in the fourth Galactic quadrant, but towards the associa- 
tion in the first quadrant. 

Recently, observational evidences of an annular volume of 
dense neutral matter, that supposedly would have been formed dur- 
ing the collision of the two bubbles, was found on the ROSAT all- 
sky surve y data byl Egger & Aschenbach' f 1995). Remarkably, as 
noticed bv lPame et al. (^1982), the molecular gas from / ~ 360° to 
290° and —25° 6 ^ 25° seemingly define a large complex of 
clouds around 100-150 pc from the Sun, comprising, e.g., p Oph, R 
CrA, G317-4, Musca, Coalsack, Chamaeleon and the Lupus clouds. 

In this case a schematic drawing of the interaction zone be- 
tween the two bubbles and the dark clouds towards the Sco-Cen 
association would be like Fi g.H] The pos ition and sizes of the dark 
clouds were obtained from Dame et al. (1987) and the ring-like 
contours from Eager & Aschenbach ( 1995). The dashed square de- 
lineates the region surveyed in this work. 

To investigate the kinematics of the interstellar gas compo- 
nents towards the Coalsack, Chamaeleon and Musca dark clouds 
we have obtained high-resolution (R ~ 60 000) spectra of the in- 
terstellar Nal D absorption lines. From the radial velocities, col- 
umn densities, velocity dispersions and the known distance of the 
stars we aim to study the physical association proposed for these 
clouds. 



2 OBSERVATIONS 
2.1 The Program Stars 

In order to investigate the interstellar gas components towards the 
Coalsack, Chamaeleon and Musca dark clouds 60 B-type stars cov- 
ering the area defined by the Galactic coordinates: 308° ^ I J? 




Galactic Longituide I (°) 



Figure 1. Schematic representation of the interaction zone between the Lo- 
cal and Loop I Bubbles with the dark clouds towards the Sco-Cen associa- 
tion. T he position and size s of the dark clouds were obtained from the sur- 
vey bv lPame et alj {l98j). The ring-like contours, taken from the ROSAT 
all-sky sui-vey data bvl Egger & Aschenbach| il99 ^. represent the annulai' 
volume of dense neutral matter supposedly formed during the collision of 
the two bubbles. The dashed square delineates the region surveyed in this 
work 

294° and —22° ^ b ^ 5° h ave been selected from the photometric 
uvbyP sample bv lCorradTfe Franco ( 1995). All targets are located 
within 500 pc from the Sun and have well determined colour ex- 
cesses E{b — y) and photometric distances ( Paper L) . SAO 257142, 
SAO 258697 and PPM 377063 were observed to complement the 
data at more negative latitudes. 

The basic information on the 65 stars is listed in Table 1 . For 
each one the visual magnitude V, the Galactic longitude and lat- 
itude, the colour excess and the distance are given in successive 
columns. Based on the uvbyP data the mean errors in the colour 
excess and distance values are 0'!"018 and 15-30%, respectively. 
The positions of the stars within the surveyed area are given in 
Fig-Hi The clouds' contours defined by the thic k lines are the lowest 
opacity level of the photographic catalogue by'Feitzineer & Stiiw^ 
( 1984) and the thin one is the outer 2 K km s^^ velocity inte grated 
CO e mission contour for the Southern Coalsack jNvman et alJ 
ll98St) . 

The data sample consists of three general divisions based upon 
the star's location within the area: inside the clouds' contours, in 
their outskirts and towards the sheet-like structure. They are, re- 
spectively, indicated by the (*), (A) and (□) signs. The horizontal 
and vertical dashed lines delineate eight sub-areas, identified by the 
roman numbers, that will be used to discuss the details of the ve- 
locity structure. Note that 5 stars are within 100 pc, while 7 other 
are between 100 and 150 pc and another 7 between 150 and 200 pc 
from the Sun. 

2.2 The Instrumental Setup 

High resolution {R ^ 60000) spectra of the Nal Di and D2 ab- 
sorption lines superposed on the continua spectra of the selected 
stars have been obtained with the 1 .4m Coude Auxiliary Telescope 
(CAT) and the Coude Echelle Spectrometer (CES) at the European 
Southern Observatory (ESO), in La Silla (Chile). The observations 
were carried out during 5 nights in April 1996 and 6 nights in April 
1997, using the remote control facilities at the ESO's HeadQuarters 
in Garching bei Munchen (Germany). 
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Table 1. Basic information on the program stars 



Star 


V 




bC) 


E{b-y) 


d(pc) 


Slai- 


V 


H°) 




E(b-y) 


d(pc) 




Towards the dark clouds contours 








Outskirts of the dark clouds 




2jyj /U 


8.754 


294.5 


3.1 


0.132 




239779 




298.1 


3.6 


0.107 




239702 


7.160 


297.6 


2.6 


0.079 


379 


240265 


7.183 


302.5 


3.2 


0.106 


489 


240621 


8.621 


305.3 


3.0 


0.274 


454 


240368 


4.617 


303.3 


3.7 


0.014 


136 


240645 


4.624 


305.5 


2.8 


0.022 


65 


251460 


7.167 


295.0 


-3.7 


0.199 


312 


251530 


8.217 


294.7 


-04 


0.105 


285 


251554 


8.215 


296.1 


-4.2 


0.272 


291 


251545 


8.018 


295.3 


-1.9 


0.115 


417 


251625 


7.332 


296.9 


-3.2 


0.149 


347 


251582 


6.880 


296.0 


-2.7 


0.164 


291 


251717 


4.724 


297.6 


-0.7 


0.072 


174 


251837 


7.519 


298.9 


l.I 


0.090 


309 


251802 


7.111 


299.2 


-3.9 


0.103 


422 


251841 


4.053 


299.3 


-1.3 


0.006 


99 


251894 


7.897 


300.6 


-6.7 


0.083 


348 


251979 


8.378 


301.3 


1.2 


0.139 


483 


251903 


4.868 


300.1 


-0.3 


0.015 


124 


251988 


8.838 


301.9 


-6.8 


0.139 


402 


251928 


6.920 


300.2 


0.8 


0.096 


313 


252114 


8.092 


303.7 


-0.7 


0.226 


498 


251987 


6.266 


301.7 


-3.6 


0.021 


117 


252245 


6.674 


305.3 


-5.7 


0.114 


123 


252262 


7.258 


305.7 


-4.6 


0.050 


172 


252284 


4.529 


306.7 


1.6 


0.018 


96 


252294 


7.810 


306.0 


-4.6 


0.108 


441 


256763 


7.978 


296.2 


-15.8 


0.465 


229 


252321 


6.198 


306.1 


-6.9 


0.108 


196 


256798 


7.675 


296.6 


-14.5 


0.220 


198 


256789 


7.970 


296.2 


-14.1 


0.161 


227 


256849 


6.883 


297.3 


-11.8 


0.286 


234 


256900 


8.113 


299.3 


-9.9 


0.194 


211 


256924 


4.240 


301.3 


-16.5 


0.007 


86 


256980 


8.440 


302.6 


-8.7 


0.159 


460 


256955 


3.841 


30 14 


-9.3 


0.003 


97 


256953 


7.777 


301.9 


-15.3 


0.160 


284 


Towards the sheet 


239327 


8.684 


294.5 


1.6 


0.068 


500 


256828 


6.399 


295.9 


-10.3 


0.157 


196 


239423 


8.239 


294.8 


4.0 


0.036 


411 


256834 


5.562 


296.3 


-10.5 


0.192 


263 


239531 


7.455 


296.2 


2.5 


0.056 


459 


256843 


6.079 


296.6 


-10.6 


0.237 


135 


240041 


7.758 


300.5 


2.8 


0.123 


231 


257002 


9.163 


303.7 


-10.6 


0.148 


368 


251629 


8.121 


296.2 


-0.0 


0.105 


433 


257032 


7.628 


305.5 


-7.8 


0.093 


284 


251699 


5.896 


298.5 


-6.7 


0.006 


99 


257033 


6.045 


305.3 


-94 


0.145 


140 


251774 


9.160 


298.1 


0.2 


0.107 


384 


257047 


6.618 


305.3 


-11.9 


0.020 


175 


251874 


8.124 


300.2 


-5.0 


0.157 


322 


257068 


8.747 


306.3 


-10.2 


0.126 


276 


251942 


8.072 


300.5 


1.8 


0.094 


454 


257069 


6.308 


305.9 


-13.0 


0.069 


149 


252136 


9.000 


304.1 


-3.2 


0.078 


353 


257142 


5.000 


307.0 


-18.0 




169" 


256745 


5.964 


297.1 


-18.3 


0.076 


255 


258697 


6.800 


306.0 


-20.3 




270" 


256788 


6.759 


294.7 


-11.1 


0143 


190 


377063'' 


5.100 


295.0 


-21.0 




165" 


256800 


5.585 


294.5 


-9.7 


0.106 


243 















" Ba.sed on trigonometric parallax from Hipparcos (ESA 19971 
^ Identification from the PPM catalogue. 



The CES optics was optimized to give a reciprocal dispersion 
of 1.8 A/mm at the central wavelength of 5890 A, resulting in a 
spectral coverage of nearly 60 A. The entrance slit width of 459 
^m, used for all the observations, yielded an actual instrumental 
resolution (FWHM) AA = 105 mA or a velocity resolution Ad 
= 5.3 km s^^ (FWHM). This corresponds to an actual resolving 
power of R ~ 56000, as determined by measurements of the in- 
trinsic width of the thorium lines. 

The CCD was a LORAL/LESSER 2688 x 512, with 16 pre- 
scan pixels. It has a pixel size of 15 x 15 /im and is relatively 
free of blemishes. Since the optical train selects only a single order, 
it was always possible to position the spectrum on a clean portion 
of the CCD. The dark current and the readout noise are very low, 
1.8 e^/pixel/hour at 164 K and 8.3 e~, respectively. The combined 
efficiency of the CAT/CES using the long camera is about 8% in 
the used central wavelength. Additional information regarding the 
CAT/CES can be found in Kaoer & Paauini ( 1996). 

Multiple sets of bias, dark and flat-field exposures were ob- 
tained each night for calibration purposes. Exposures of a thorium 
lamp were taken before and after each object exposure to ensure a 
better wavelength calibration. The CES, however, showed to be a 
very stable instrument. 

Exposure times ranged from 2 minutes to 2.5 hours, depending 
on the visual magnitude and spectral type. In all cases signal-to- 



noise ratios (S/N) greater than 150 in the raw spectra were achieved. 
At least two exposures were obtained at each setup position and 
the integration times were limited to a maximum of 30 minutes to 
minimize the effects of cosmic ray events in the individual spectra. 



2.3 Stellar Distances 

Although Hipparcos satellite was laun ched in 1989 August, the cat- 
alogue only became available in 1997 <ESAll997h . For this reason, 
the observing sample was selected based on our photometric dis- 
tance determination. However, the paralactic distances may now 
serve as a basis for comparison with our photometric distances. A 
search in the Hipparcos catalogue provided us 5 1 stars in common 
with our sample - one of them, SAO 239370, appears with neg- 
ative value for its trigonometric parallax, and was excluded from 
the comparison. Fig.|3|shows a plot containing the distances based 
on the Hipparcos parallax vs. photometric distances. The denoted 
errors for the photometric distances were estimated for each star 
based on the individual uncertainties of their photometric measure- 
ments I Corradi & Franco 1995), and ranges from 1 1 to 49 per cent. 
In general, the agreemment is very good testifying the high qual- 
ity of our photometric distance determination. There are few stars 
showing rather large discrepancy between their estimated photo- 
metric and parallatic distances, usually those beyond 200 pc and 
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B-type stars - < d < 500 pc 




308 306 304 302 300 298 296 294 
Galactic Longitude / (") 



Figure 2. Location of the stars within the studied area. The clouds' contours 
defined by the thick lines are the lowest opacity level of the photographic 
catalogue by Feitzinger & Stiiwe 1 1984) and the thin one is the outer 2 K 
km velocity integrated CO emission contour for the Coalsack. The (*), 
(A) and (□) signs indicate the stars with line-of-sight inside the clouds' 
contours, in their outskirts and towards the sheet-like structure, respectively. 
The horizontal and vertical dashed lines dehneate eight sub-areas, identified 
by the roman numbers, that are used to discuss the details of the velocity 
structure. 

large uncertainties in their trigonometric parallax. However, as will 
be seen further, the main interstellar structures discussed in this 
work are located nearer than that, and consequently these uncer- 
tainties in the distance estimation do not affect substantially our 
conclusions. 



3 DATA REDUCTION 

The initial processing of the CCD frames employed IRAF routines 
to subtract the bias, divide by a normalized flat-field, and remove 
cosmic rays from the rows occupied by the stellar spectrum and the 
background regions. The one-dimensional spectra were then ex- 
tracted using IRAF routines. 

It is worth reiterating that usually the dome flat-fields pro- 
vided the best correction for the strong vignetting present in the 
spectra borders, besides that the affected parts were trimmed off 
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Figure 3. Comparison between the trigonometric distances obtained from 
the Hipparcos catalogue lESA 1997.) and the photometric distances based 
on uvbyP data. 

the images. The wavelength calibration was established with a set 
of 20 to 30 lines , iden tified using the thorium line wavelengths of 
iD'Odorico et alJ ^1987). A second-order Chebyshev function was 
fitted to the line wavelength as a function of the pixel number, with 
rms scatter of the residuals about the fit smaller than 0.002 A (0.05 
km s~^) in all cases. The dispersion solution was applied to all sci- 
ence exposures with no rebinning of the data, and the final spectra 
normalized to unity continuum with cubic splines. 

The individual normalized spectra were then co-added on a 
pixel-by-pixel basis, weighted by the inverse of the rms deviation 
of the continuum fit. No noticeable degradation in spectral resolu- 
tion resulted from the co-addittion procedure, since the individual 
spectra have been obtained within an hour or two of each other, and 
the velocity changes in the observer rest frame, over this time scale, 
are very small compared to the instrumental resolution. 

Numerous absorption lines arising from the Earth's atmo- 
sphere occur in this region of the spectrum, and two lines in partic- 
ular (at 5889.637 A and 5890.09 A) are very close to the rest wave- 
length of the D2 line. Hobbs r l978) discusses the origin and iden- 
tification of these atmospheric (telluric) lines in the Na I D wave- 
length region. To reduce telluric line contamination the stellar spec- 
tra were divided out by a purely atmospheric absorption template 
spectrum obtained from a fast rotating, unreddened early-type star 
with no interstellar sodium, and taken at similar air mass. For both 
observing runs spectra of a Leo and 7 Lup were taken each night. 
The equivalent width of the Nal D lines for these stars are esti- 
mated to be ^ 1 mA and 1.4 mA, respectively (Welsh et al. 199 J 
ICrawforJl99lh . The stars SAO 251050 and SAO 224833 (tLup), 
initially chosen as candidates for telluric line correction, showed 
the presence of interstellar Na I D lines. 

Since most of the telluric lines arise from water molecules, 
their strengths change throughout the night as the water content of 
the Earth's atmosphere changes. Hence, the template spectra were 
properly scaled to match the strengths of the atmospheric lines in 
the object spectra. The success of the atmospheric correction pro- 
cedure is apparent from the fact that all the telluric lines, which are 
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clear of the interstellar lines, have been satisfactorily removed (see 
Fig.EJ. 

Finally the observed wavelengths were brought to the velocity 
scale and converted to the Local Standard of Rest (LSR) velocity 
frame, a ssuming a solar motion of 20 km s~^ towards 0=18*^ and 
6 = 30° <Mihalas & Binnevll98lh . To access the zero point of the 
velocities the wavelengths of some of the numerous telluric lines 
gathered in the observed spectral range were measured, providing 
an external accuracy of the order of 4 mA and showing that the 
lines match quite well, seemingly with no systematic offset. 



4 THE ABSORPTION-LINE PROFILES 

All observed absorption lines show appreciable velocity structure, 
and thus profile fitting provides the best method for determining ac- 
curate interstellar column densities as well as to discern and deter- 
mine the properties of the individual interstellar clouds contribut- 
ing to the line profile. The observed profiles of the D2 and Di 
absorption-lines are shown in Fig. |4| The dots represent the ob- 
ject spectra and the solid line gives the adopted multicomponent 
fit, as described below. The line-fitting program (FITS6P) has been 
kindly provided to us by Dr. Welty. Fur t her de tails of the employed 
metho d can be found in IWeltv et'Sl J 1994^ and 'Sem bach et alj 
il993h . The zero point of the velocities in the upper and lower pan- 
els are with respect to the rest wavelength of the D2 (5889.9510 
A) and Di (5895.9242 A) lines, respectively. In each panel the tick 
marks and numbers above the profiles identifies the different com- 
ponents used to fit the D2 and Di lines. 

Voigt line profiles parameterized by their velocity dispersion 
parameter b, cloud component velocity Visr and cloud column 
density A^NaI> were convolved with a Gaussian instrumental profile 
of 5.3 km width and fitted to the observed data. Inspection of 
the calibration lamp spectra supports our assumption of a Gaussian 
instrumental profile since no asymmetric tails are observed in the 
thorium lines. Both the D2 and Di lines were fitted simultaneously, 
with the fewest components that seemed reasonable. Inflections in 
line shape, as a result of additional components in the wings of the 
stronger lines, often indicated whether another component should 
be added to the model. The hyperfine structure splitting of ~ 1 
kms~^ of the individual D lines were not included in the adjust- 
ments since the resolution of the data is sufficiently broad to be 
affected by the splitting. 

In most cases the assigned initial values of b were from 2.0 
to 4.0 kms~^ to all components, but larger values are sometimes 
required, e.g. SAO 239531. Also the three free parameters (A^NaI> 
Visr , b) for each component were allowed to vary in an iterative 
nonlinear least-squares fit to the observed profile. Some parame- 
ters occasionally were held fixed at "reasonable" values in order 
to facilitate convergence of the fits, but are thought to be as well 
determined as those allowed to vary. 

The average of the D2 and Di best fit values of these param- 
eters are listed in Table 2. A total of 125 entries are given. Each 
successive column gives the star's name, component number, LSR 
radial velocity Viar , heliocentric radial velocity Vq, logarithm of 
the column density log N-f^^ i , velocity dispersion b and the equiv- 
alent widths Wd2 and Wdi, respectively. 

The distribution of the absolute mean errors for log A'^^a I > 
Visr and b are shown in Figs.|5^-c. For narrow, isolated lines of 
moderate strength such errors are typically 1%-10% for log A'^Nal 
and few hundredths of a km s~^ for b and Visr ■ For stronger, wider 
and more blended lines the uncertainties in logA'^Nal srs typi- 



cally 5%-20%, and 0. 1-0.4 km s~^ in b and Visr ■ For the strongest 
saturated lines, uncertainties in the actual component structure are 
likely to be the dominant error source. 

The theoretical uncertainty in the placement of the continuum 
level is approximately equal to the rms scatter divided by the square 
root of the number of points used to define it ( Howarth & Phillipa 
1986). The statistical error on line flux values are given by ap = 
^s/n" ' where n is the number of observed data points in the ab- 
sorption profile, A A is the resolution element in A/pixe l, and S/N 
is the signal-to-noise ratio of the exposure (Welsh et alJ I 199(1) . In 
most cases the continua are very well defined and statistical errors 
probably dominate the equivalent width errors. Typically for our 
observations, n = 16, A A = 0.027 A/pix and S/N ^ 150, so that 
cTp ^ 1 mA. 

The star tLup has been previo u sly ob served at similar res- 
olution (4 kms~^) by IWelsh et alj il994h . The parameters for 
(logiVNal, Visr,h) are (10.9 ± 0.02, -13.8 ± 0.3, 4.6 ± 0.3), 
and those obtained in this work are (11.0 ± 0.3, -13.42 ± 0.09, 
4.3±0.7). The remarkable aggreement between the values is rather 
encouraging. Additional absorption features occur in some of the 
spectra of low v sin i stars shown in Fig.|4| The stellar Na I D lines 
are expected to be very weak throghout the class of B stars (being 
^ 2 mA for stars earlier than B7), but the cases where stellar lines 
have been identified were removed from the multicomponent fit. 



5 THE KINEMATICAL STRUCTURE 

Histograms of the individual distributions of log N-^^^ i , Vjsr and b 
are given in Figs.|5jl to|5f. Inspection of Fig.|5li shows a clear divi- 
sion of the column densities around log A'^^al ~ 12.0 cm^^ . The 
average column density of the sub-samples is 11.2 cm~^ and 12.7 
cm^'^ , respectively, both with small dispersion of ± 0.6 cm~^ . 
Figure [5^ shows that the distribution of the LSR velocities, al- 
though centered around kms~^, is skewed to the negative side. 
The velocity dispersion parameter b shows a concentration around 
2.5 kms^^ and barely suggests another around 5.5 kms^^. These 
trends are somewhat confirmed by Figs.|5^-i, that shows the plots 
of b tis. logA^Nal. logA^Nal 'VS. Visr and Visr fs. b, respec- 
tively. For the sake of clarity few points of very high negative and 
positive LSR velocities have been excluded from these plots. 

As one can see the lower column densities span the whole ve- 
locity dispersion range, and are concentrated mostly on the negative 
velocity range. On the other hand, the higher column densities are 
concentrated around b ~ 3 kms~^, and have velocities spread by 
± 3 kms~^ around zero kms~^ in the LSR frame. The apparent 
lack of points with log A^^al ^ 10-5 cm~^ for all b may be due 
to our detection limits. The lack of points with b ^5 1.5 kms^^ 
for log Nj^^ I 11 .5 cm^'^ may be due to our inability to resolve 
closely blended components. 

5.1 Distribution of the interstellar gas components 

To investigate the distribution and kinematics of the interstellar gas 
components along the line-of-sight plots of logTV^ai, Visr and 
b as a function of the stellar distance are given in Fig.|6l The left 
column shows those stars with line-of-sight inside the clouds' con- 
tours; the middle column those in their outskirts; and the right col- 
umn those towards the sheet. 

The plotting simbols indicate a further division of the compo- 
nents according to the line profile characteristics. The line profiles 
with a single weak feature are represented by the open triangles. 
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Figure 4. The observed Na I D absorption Unes. The observed intensities are plotted as dots; the smooth curves are theoretical line profiles with the parameters 
given in Table 3. The tick marks and the corresponding number identifies the fitted components. The zero point of the velocity scale is referred to the rest 
wavelength of the D2 and Di lines 



while those profiles with a single strong feature are represented by 
the solid triangles. The profiles with two weak features are repre- 
sented by the plusses and asterisks, which indicate the features with 
velocity similar to the single weak and single strong components, 
respectively. The profiles with two or more features are represented 
by the squares and circles. The open and solid squares indicate 
those features similar to the single weak and to the single strong 
components. The open circles indicate the other components, i.e., 
the third, fourth and so on. 

The structure of the colunm density vs. distance diagrams for 
the stars inside the clouds' contours, in their outskirts and towards 



the sheet is fairly similar. Single weak components are seen for the 

stars closer than 60-80 pc from Sun (open triangles and plusses), 
and for most of the more distant stars along the distance range 
(open squares). On the other hand, the single strong component 
is seen only for stars more distant than 120-150 pc from the Sun 
(solid triangles, asterisks and solid squares). The open circles tend 
to appear at distances greater than 300-350 pc, but without clear 
correlation with the column densities or the radial velocities. 

The velocities of the components identified along the line of 
sight also have a similar structure for the stars inside the clouds' 
contours, in their outskirts and towards the sheet. The single weak 
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Figure 4. (Continued) 



component (open triangles, open squares and plusses) has velocity 
more negative than —5 km s~^, while the single strong component 
(solid triangle, solid square and asterisks) is centered at km s^^ in 
the LSR frame. The fact that the plusses and asterisks show kine- 
matic structure consistent with those of the single weak and the 
single strong components, respectively, suggests that the material 
have some inhomogeneities. The open circles appear in the whole 
range of radial velocities. 

The velocity dispersions of the components also have a similar 
structure for the stars inside the clouds' contours, in their outskirts 
and towards the sheet. The single weak component has b ~ 5.5 
kms~^. On the other side the single strong component has b ~ 2.5 
kms~^. 



The b value is a measure of both thermal and turbulen t veloc- 
ities in an interstellar cloud. We have b — a/^^ + 2v'}, where 
vt is the rms turbulent velocity, k is the Boltzmann's constant, T 
is the kinetic temperature and m is the mass of the atom. There is 
no independent way to determine the degree of turbulence within 
the material, but some qualitative considerations can be done. The 
typical kinetic temperature of interstellar diffuse clouds is gener- 
ally taken to be a bout 80 - 115 K, fr om observations of rotational 
excitation of Ho <Savage et al1ll97'1) . For Na I this corresponds to 
0.24 ^ b ^ 0.29 kms~^. In this case the upper limits to vt for 
the weak and strong components would be 1.8 and 3.9 kms^^. 
Assumi ng a sound sp eed for diffuse interstellar clouds about 0.7 
kms" jSDitzeJl978l) the observed components are either signifi- 
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Figure 4. (Continued) 



cantly hotter than a standard diffuse cloud or are subject to super- 
sonic turbulent motions. 



5.2 The gas sheet-like features 

To loolc for a possible association of the interstellar gas components 
with the observed dust features in this direction, log A^^j^j , Visr 
and b as a function of the E{b — y) colour excess have been plotted 
in Fig.|7| The columns and symbols have the same meaning as for 
Fig.|6| It can be seen that the strong component is only picked up 
by the reddened stars, i.e., those with E(b — y) 0™050, whereas 
the weak component is picked up by both reddened and unreddened 



stars, suggestting that dust and gas apparently have the same distri- 
bution along the line of sight. 

The average column density of the weak component 
(log A'^Nal ~ 11.2 cm"^ ), would correspond to iVn ~ 3 x 10^^ 
cm~^ . It is worth noticing at this point that, in Fig. 6a of Paper 
I, around 60-70 pc the dispersion of the E{h ~ y) increases, with 
the lower limit becoming slightly shifted to positive values. This 
indicates that the spectroscopy has provided a finer tuning of the 
existing low column density components, more difficult to disen- 
tagle based only on the colour excess values. On the other hand, it 
shows that the observed increase in the dispersion of E[h — y) is 
real, and not due to observational errors, attesting the high quality 
of the photometric data and method. 
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Figure 4. (Continued) 



There appears to be a common envelope of minimum and 
maximum column density, centered in a narrow distance slot 
around 120-150 pc, that remain unchanged up to 350-400 pc, 
suggestting that a second volume of low density may have been 
reached. The complete range of column densities from log A^^a I ~ 
10.7 to 1 1.7, is also shown by the stars centered around 60 pc, sug- 
gesting a sheet-like structure. 

Interestingly, the minimum and maximum column densities 
of the strong component (log A'^^al = 12-3 and 13.2 cm^^ ) cor- 
respond to E{b - J/) ~ 0'"050 and 0'"300, respectively, if it is 
assumed that the standard relationship between log N-f^^^ j and log 
A^'h (perlet et al. 1985) and E{b - y) and log Nh (Knude 1978) are 
valid. This may suggest that the high column densities may be the 



gas related to the dust sheet-lilce structure, reported in Paper I. Since 
the Coalsack, Musca and Chamaeleon are also located around 150 
± 30 pc they seem also to be part of the gas composing the sheet. 

Assuming that the observed stars are within different dark 
clouds the column densities, radial velocities and velocity disper- 
sions, shown in Figs.|4|to|6| could be either the effect of a common 
structure containing the embedded molecular clouds or just an ef- 
fect of different superposing structures. In order to investigate how 
the parameters ( log Nj^i^ j , Visr ,b) are distributed over the studied 
area we have used the eight sub-areas delineated by the horizontal 
and vertical dashed lines and identified by the roman numbers in 
Fig.H 

Figures [8^ and b show the distribution of logA^^j^i (left). 
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Figure 4. (Continued) 



Visr (middle) and b (right) against the Galactic longitude and lat- 
itude, respectively, for the eight sub-areas. The symbols have the 
same meaning as for Fig.|6| and the vertical and horizontal dashed 
lines indicate the division between the sub-areas. The strong and 
weak components are present - and well separated - everywhere. 
The gas components form two extended interstellar features, that 
cover the whole connecting area. 

Note, however, that the radial velocity of the strong compo- 
nent seems to increase from —3 kms^^ near the galactic plane 
to -1-3 kms~^ near the southern parts. This trend would be ex- 
pected from an expanding bubble with the center around I ~ 309°, 
b — —9°. This fact is quite interesting given the fact that high ther- 
mal pressure has been measured for the hot gas phase for a nearby 



direction (I = 309°, b — —15°) from a combination of ROSAT 
0.25 keV and optical data (Knude et al. - in preparation). 



6 DISCUSSION 

The analysis of the identified components indicates that the inter- 
stellar gas is distributed in two extended sheet-like structures per- 
meatting the whole area, one at d ^60 pc and another around 
150 pc from the Sun. The nearby feature is approaching the Sun 
with average radial velocity of —7 kms~^, has low average col- 
umn density ( log A'^^a I ~ 11-2 cm~^ ) and is either hotter or more 
turbulent than the strong component, with velocity dispersion b ~ 
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Figure 4. (Continued) 



5 kms~^. The more distant feature lias column densities between 
12.3 ^ logA'^Nal ^ 13-2, average velocity dispersion b « 2.5 
kms^^ and seems associated to the dust sheet observed towards 
the Coalsack, Musca and Chamaeleon direction. Its velocity is cen- 
tered around kms~^, but there is a trend for increasing from —3 
kms~^ near 6 = 1° to 3 kms~^ near b = —18°. 

There seems to be a common envelope of minimum and max- 
imum column density, centered in a narrow distance slot around 
120 - 150 pc, that remains unchanged up to 350 pc, suggestting 
that a volume of low density may have been reached. After that, 
a third group of high negative and positive velocity components 
seemingly arise around 300-350 pc from the Sun, however, without 
indications of an extended structure. 



6.1 The nearby, low column density feature 

Details of the overall distribution of the ISM (d ^ 500 pc) along the 
studied direction have been extensively discussed i n Paper I. The 
kinematics of the nearby gas has been reviewed byl prisch & York! 
( 1986). Details of the velocity structure in the region (360°^ I ^ 
295° and 0°^ b ^ 30°), encom passing most of the Sco-Cen asso- 
ciation, have been discussed bv iCrawford il99il) . 

In the immediate solar neighbourhood the observational data 
suggest that the Sun is located within but close to the edge of 
the Local Interstellar Cloud (LIC). First identified from its kine- 
matics by |L allem ent & Berlin ( 1992), the LIC was modelled by 
iRedfield & Linskvl [200Cl) as roughly spherical with dimensions of 
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Figure 4. (Continued) 



5-8 pc, warm (T fa 7000 K), low density (nni < 0.1 cm~^) and 
mostly neutral (ric/nm < 0.5). Its largest column density is about 
2 X 10^** cm"^ towards / = 157°, b = -25°. 

Moving at a helioce ntric velocity of ~ 26 kms~^ toward s 
I « 186°, b f» -16° (cf. Lallement et al. 'l 995"; 'Lallement" 1998"), 
the Lie flows through the solar system and resonantly scatters solar 
Hi (Lyman a) and Hel (A584) radiation. This flow, is referred to as 
the local interstellar wind (LISW). 

The Lie and other warm clouds with similar temperatures but 
a wide range of metal depletions, are embedded in an irregularly 
shaped region, whose radius ranges from 30 to 300 pc, and is de- 
ficient i n dense neutral hydrogen compared to the galactic aver- 
age fe.g. lFrisch & Yorldll983t[Parescdll984ISnowden et alll990l: 



lTinbergerj[l98S IWarwick et al]|l993t IWelsh etalJll994 . Usually 
called Local Bubble, this cavity is thought to contain hot (T ~ 
10*" K), low-density (tihi ^ 0.025 cm^'^) gas and as recently 
shown by the ROSAT and EUVE shadowing experiments, coex- 
ists with the neutral atomic and molecular gas within its interior 
(e.g. feowver et al. 1995; Kerp et al. 1993; Snowden et al. 199ll 
1995, 1998; Wane & Yu 1995). It is, however, worth noticing that 
iCravens 12000.) has recently shown that heliospheric X-ray emis- 
sion can account for about 25%-50% of the observed soft X-ray 
background intensities, which may put in doubt if the Local Bub- 
ble really contains the million degree hot gas. 

Optical interstellar absorption line studies within 100 pc of the 
Sun show evidence of a flow of material analogous to the LISW. 
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Table 2. Profile fit results for the Na I D absorption-lines. Each successive column gives the star's name, component number, LSR radial velocity Vi^^ , column 
density ( log A^n^ i ), velocity dispersion b and the equivalent width of the D2 and Di lines 
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Assuming that the Sun is embedded within a single coherently 
moving in terstellar wind, from a direction (Zw,bw) with velocity Ww, 
ICrutcheilh982 ^ obtained the following vect or (Z„, h^,v^) = (345° , 
-10°, -15 kms~^). In a similar anal ysis iFrisch & Yorld <1986t) 
obtained (354°, +3°, -12 kms'^) and lCrawford Jl99lh obtained 
(301°, +59°, —9 kms^^). A possible description for this flow as- 
sumes that the LISM is material which belongs to the surface of 
an expanding shell centered on the Sco-Cen association. ICrawford 



il99ll) has obtained an expansion velocity of 7 - 9 km s^^ assum- 
ing the center at (I = 320°, b = +10°), the dist ance to the Sun of 
140 pc and the radius of 1 10 pc, as suggested bv ldeGeujfl993) . 

Another possible explanation for the existence of the neutral 
clouds and their predominant flow away from the Sco-Cen asso- 
ciation is a consequence of the interaction between the Local and 
Loop I Bubbles (see Se ct.l6.2> and subsequent loc al fragmentation 
of the interaction zone iSreitschwerdt et all200(l) . 
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Figure 5. Distribution of the free pai'ameters determined in the profile fits and their absolute mean eiTors. a-c Distiibution of the absolute mean errors in the 
fitted values of logA'^njal. Vlar and b, respectively. For narrow, isolated lines of moderate strength such errors are typically 1%-10% for logA^Nai and few 
hundredths of a kms~^forband Vj^^ . For stronger, wider and more blended lines the uncertainties in logA^N^I are typically 5%-20%, and 0.1-0.4 kms~^ 
in b and Visr ■ d-f Histogram of log A^n^ i , Vj^^ and b, respectively. There is a clear division of the column densities around log A^n^ i ~ 12.0 cm~^ . The 
LSR velocities, although centered around km s~^, are skewed to the negative side. Also the velocity dispersion parameter b shows a concentration around 
2.5 kms~^ and another around 5.5 kms~^. (g) Distribution of logA^j^jjii vs. b. (h) Distribution of log A'^f.jai i'.v. Vja^ . (i) Distribution of Vj^r vs. b. The 
lower column densities span the whole velocity dispersion range, and are concentrated mostly on the negative velocity range. On the other hand, the higher 
column densities are concentrated around b « 3 km s^^, and are essentialy at rest in the LSR frame; although a ± 2 km s^^ spread may be present. For the 
sake of clarity few points of very high negative and positive LSR velocities have been excluded from these plots. Further details can be found in the text 



Since the closest star in our sample is at 65 pc the LIC can- 
not be addressed properly, neither the distance of the lower column 
density component. However, according to Welsh et al. ( 1994, and 
references therein) the column densities of the nearby feature ob- 
served in this work should not be much closer than 50 pc from the 
Sun. 

For the latitude/longitude region observed in our work the ve- 
locities predicted by the interstellar wind flow vector a nd the ex- 
pandi ng shell models give essentially the same results <Crawford 
I1991L cf. his Fig. 3). In any case the velocity of the low column 
density component is consistent with the previous results, and indi- 
cate an outflow from the Sco-Cen association. 

A more recent discussion bv lChaetalJ j200dl) also using Na I 
D line data and Hipparcos distances proposes components at three 
distinct heliocentric velocities in the Puppis-Vela region. Their 



component A with the velocity in the range from -1-6 to -1-9 km s 
covers part of the region presently discussed. 

From the distances of Table 1 and heliocentric velocities from 
Table 2 w e construct Fig. |9| of Vq vs. distance and notice that 
ICha et alT s component A is seen in two stars closer than 100 pc and 
in fact seems present in all stars all the way to 450 pc, this implies 
that component A covers all of our surveyed region, indeed a sheet 
nearer than about 100 pc. Two stars marginally have absorption fea- 
tures at component B's lower velocity limit (12 km s^^), these stars 
are at 100 pc and 130pc respectively. Otherwise component B is 
seen in most stars beyond 180pc; the upper distance limit is very 
sharply defined by component B's upper limit at 15 kms^^. No 
stars beyond ~350pc shows any indication of component B. Com- 
ponent C (21 - 23 kms~^) is virtually absent from our re gion an d 
is only observed for galactic longitude less than 275° in Ch a et alf s 
survey, far from our region. Four stars closer than 100 pc show ab- 
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Figure 6. Plots of log A^n^i , Vi^r and b against the stellar distances. From left to right the columns show the stars with line-of-sight inside the clouds' 
contours, in their outskirts, and towards the sheet. The line profiles with a single weak feature are represented by the open triangles, while the profiles with a 
single strong feature are represented by the solid triangles. The profiles with two weak features are represented by the plusses and crosses, which indicate the 
features with velocity similar to the single weak and strong components, respectively. The profiles with two or more features are represented by the squares 
and circles. The open and solid squares indicate those features similar to the single weak and strong components. The open circles indicate the other ones, 
i.e., the third, fourth, etc. The saine structure for the weak and the strong components suggests one interstellar gas at d ^ 60 pc and other 120-150 pc. The 
strong coinponent (logA^N^i fa 12.7 cm~^ ) is spread ± 3 kms~^around the zero velocities, and has low velocity dispersion b Ri 2.5 kms~^. The weak 



component ( log A^Na I ~ 
Further details in the text 



11.2 cm ) is approaching to the Sun with average Vig^ around —7 kins and have larger velocity dispersion b = 5 kins 



sorption in a narrow interval from +1 to +4 km s~^, this interval is 
hardly represented between 100 and 400 pc from where it is again 
present in almost any star. A fifth group of absorption has veloc- 
ities between —2 and —4 kms^^ for three stars between 120 and 
140 pc, lines in this interval disappear for larger distances but reap- 
pear between 300 and 400 pc. Locally, that is for distances less than 
150 pc, we confirm three interstellar kinematic groups with absorp- 
tion lines in the intervals: [—5,0], [0,5] and [5,10] the second of 
which represents the volume within lOOpc whereas the other two 
both are present between 120 and 1 50 pc. The local ISM thus seems 
more complicated than proposed in IChaetalJ 12000) investigation. 
One could speculate whether the interval [0,5] is representing the 



local undisturbed material and [5,10] and [—5,0] the expansion of 
the Local and the Sco-Cen Bubbles respectively? 

Fig- m also indicates that the distribution of heliocentric ve- 
locities for the diffuse components changes character between 300 
and 350 pc. Closer than 300 pc there is a set of velocities exceeding 
10 kms~^, this set disappears beyond 350 pc. On the other hand a 
new set with heliocentric velocities smaller than km s~^ appears 
beyond 300 pc. Could this kinematic transition at 300-350 pc be 
connected to the expansive movement of the backside of the Lower 
Centaurus-Crux feature being countered by an opposite motion? 
The expansion has ceased at about 350 pc and the counter motion 
halting the expansion is first noticed at about 300 pc. 
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Figure 7. Plots of log Afj^^I > Vlsr t> against the E{b — y) colour excess. From left to right the columns show the stars with line-of-sight inside the 
clouds' contours, in their outskirts, and towards the sheet. The plotting symbols have the same meaning as in Fig.|6| The strong component is only piclced up 
by the reddened stars, whereas the weak component is picked up by both reddened and unreddened stars. Such diagrams suggest that dust and gas apparently 
have the same distribution along the line of sight 



6.2 The interface between the Local and Loop I bubbles 

A soft X-ray shadow has recently been discovered on the edges 
of the Loop I Bubble, cast by a warped annular volume of dense 
neutral matter that supposedly forme d duri ng its collision with the 
Local Bubble llEgger & Aschenbachlll995i) . The shadow counter- 
part, a huge HI ring was also identified on the lPickev & LockmanI 
il99d) data. 

lEgger & AschenbachI il995h suggested that the steep increase 
of the column density caused by the annular feature, from less 
than 10^° to ^ 7 xlO'^'^, occurs at a distance of ~ 70 pc from 
the Sun. This distance is thought to be supported by results from 
optical and UV spectral analysis of stars nea r the center of Loop I 
(310° ^ I ^ 330° and 15° sC 6 ^ 25°) bv lCenturion & Vladilol 
il99lll . who found the presence of a neutral gas wall of A^h ~ W^" 
cm~^ at a distance of 40 ± 25 pc. See also the determination of 
the radial extent of the Local Hot Bubble recently estimated by 



ISnowden etall h998h based on the ROSAT and IRAS all-sky sur- 
vey data. 

However, for the hydrogen column density of the annular fea- 
ture, the colour excess and the Nal column densities obtained in 
this work suggest that the interaction zone between the two bub- 
bles is locate d around 120-150 pc from the Sun. This result is cor- 
roborated by 'Sfeiretal.71999') - an improved version of the 3D 
mapping of the L ocal Bubble cavity has been recently obtained 
bv lLallement et alJ l2003l) - where it is found that for the region 
investigated here the 20-50 mA isocontours at b ~ —20° occurs 
at d ~ 100 — 150 pc, but at 6 ~ —45° the distance is close to 
70 pc. The hydrogen column densities and velocities observed by 
ICenturion & Vladilo (1991) are more consistent with our lower col- 
umn density component at d ^ 60 pc. In the same strip of latitude, 
the wall of neutral gas has a column density A'^h ~ 10^^ cm~^ , 
as suggested by Iwan ( 1980). Moreover, the hydrogen column den- 
sities towards the interaction ring, used bv lEgger & AschenbachI 
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Figure 8. a Distribution of log A^nj^i (left), Vi^r (middle) and b (right) with the Galactic Longitude. The symbols have the same meaning as for Fig. 1^ and 
the dashed lines indicate the dividing point between the sub-areas. The weak and the strong components form two extended interstellar components, that cover 
the whole connecting area. There seems to have no dependence with the Galactic longitude, but the radial velocity may be increasing with the latitude 



i 19951 cf. their Fig. 4) suggest a value around 3 x 10^^ cm ^ 
at 70 pc, also more consi stent with our nearby feature. Since 
ICenturion & Vladild iT991 ^ have observed only unreddened direc- 
tions our results cast some doubt on the distance to the interaction 
zone being located around 70 pc. It seems that the annular region 
is twisted and folded, with different directions having different dis- 
tances. 

In addition to the velocities indicative of the flow of material 
in the last section, some of the more d istant stars also have compo- 
nents at nearly kms"^ in the LSR <Frisch & York 1986). In the 
longitude range 360 ^ / ^ 345 there are a number of low velocity 
(— 4 ^5 Visr ^ -f 4 km s^ ^ ) components with generally strong Na I 
lines, consistent with the proximity of the p Oph and Lupus clouds 
(Crawford 1991). In the longitude range 325 ^ I ^ 295 there are 
also several components, in the positive latitudes, with LSR veloci- 
ties which are small but clearly negative (— 4 ^ Vjsr ^ km s~^). 

As mentioned in Sect.|5| the feature around 120 - 150 pc ob- 
served in our work, has column densities suggestting that the gas is 



associated to the dust sheet observed towards the Coalsack, Musca 
and Chamaeleon direction. It is certainly not appropriate to make 
generalizations exclusively from our data, but the fact that the ve- 
locities observed here are also in the same range (—3 ^ Visr ^ +3 
kms~^), suggests that the dust and gas feature around 120 - 150 
pc, seems to be part of an extended large scale feature of simi- 
lar kinematic properties, supposedly identified with the interface of 
the Local and Loop I bubbles. 



7 CONCLUSIONS 

The investigation of the interstellar gas components towards the 
Southern Coalsack, Chamaeleon and Musca dark have produced 
the following results: 

• The interstellar gas is distributed in two extended sheet-like 
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Figure 8. b Distribution of log N-^^ \ {left), Visr {middle) and b (right) with the Galactic Latitude. The symbols have the same meaning as for Fig.|6|and the 
horizontal dashed lines indicate the dividing point between the sub-areas. The weak and the strong components are present in the whole latitude range. Note 
that radial velocity of the strong component seems to increase from —3 kms~^ near the galactic plane to -h3 kms~^ near the southern parts 
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structures permeatting the whole searched area, one at d ^ 60 pc 
and another around 120 - 150 pc from the Sun. 

• The nearby feature is approaching to the Sun with average 
radial velocity of —7 kms~^, has low average column density 
log7Vi,jai~ 11-2 cm~^ and has velocity dispersion b ~ 5 kms~^. 
The more distant feature has column densities between 12.3 ^ 
log A''i>ja I 5; 13.2, average velocity dispersion b ~ 2.5 km s^^ and 
seems associated to the dust sheet observed towards the Coalsack, 
Musca and Chamaeleon direction. Its velocity is centered around 
kms~^, but there is a trend for increasing from —3 kms~^ near 
& = 1° to 3 kms~^ near 6 = -18°. 

• In aggreement with several independent investigations the 
nearby low column density feature indicates a general outflow from 
the Sco-Cen association, while the dust and gas feature around 120 
- 150 pc seem to be part of an extended large scale feature of simi- 
lar kinematic properties, supposedly identified with the interaction 
zone of the Local and Loop 1 bubbles. Assuming that the inter- 
face and the ring-like volume of dense neutral matter have similar 
properties, our results suggest that the interaction zone between the 
bubbles is twisted and folded. 
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Figure 9. Heliocentric velocities of the interstellar clouds components 
as a function of stellar distance. Filled circles represent components for 
which logA^Nai ^ 12.5 cm^^, while open circles represent components 
for which \ogN^^\ < 12.5 cm~^. Note the well defined velocity distribu- 
tion for the strong component and that it appears for stars beyond ~ 140 pc. 
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